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Abstract. We study the contribution of projection effects to the intrinsic thickness of the Fundamental Plane (FP) 
of elliptical galaxies. The Monte-Carlo mapping technique between model properties and observed quantities, 
introduced by Bertin, Ciotti, & Del Principe (2002), is extended to oblate, two-integrals galaxy models, with 
non-homologous density profiles, adjustable flattening, variable amount of ordered rotational support, and for 
which all the relevant projected dynamical quantities can be expressed in fully analytical way. In agreement with 
previous works, it is found that projection effects move models not exactly parallel to the edge-on FP, by an 
amount that can be as large as the observed FP thickness. The statistical contribution of projection effects to the 
FP thickness is however marginal, and the estimated physical FP rma thickness is ~ 90% of the observed one 
(when corrected for measurement errors). 
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1. Introduction 

In the observational three-dimensional space of central 
velocity dispersion oo, (circularized) effective radius R e , 
and mean surface brightness within the effective radius 
7 C , early-type galaxies approximately locate on a plane, 
called the Fundamental Plane (hereafter FP; Dressier et 
al. 1987; Djorgovski & Davis 1987), and represented by 
the best-fit relation: 



log R c = a log cr + (3 log I c + 7. 



(1) 



The coefficients a, (3, and 7 depend slightly on the con- 
sidered photometric band (e.g., Pahre, de Carvalho, & 
Djorgovski, 1998; Scodeggio et al. 1998). By measuring 
R c in kpc, (To in km/s, and I c = L/(2ttR^) in L Q /pc 2 
(where L is the total galaxy luminosity), reported values 
in the Gunn r band are a = 1.24 ±0.07, (3 = -0.82 ±0.02, 
7 = 0.182 1 (J0rgensen, Franx, & Kjaergaard 1996, here- 
after JFK96). One of the most striking observational prop- 
erties of the FP is its small and nearly constant scatter: 
the distribution of logi? e around the best-fit FP has a 
measured rms (after correction for measurement errors, 
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1 This value of 7 refers to the Coma cluster and to Ho = 50 
km s -1 Mpc" 1 . 



hereafter (Tint), that corresponds to a scatter in R at fixed 
(T and / ranging from 15% to 20% (see, e.g., Faber et al. 
1987; JFK96). 

For a stationary stellar system the scalar virial theorem 
can be written as 



Re 



= Kva, 



v^o, 



(2) 



where T* is the stellar mass-to-light ratio in the pho- 
tometric band used for the determination of L and i? e , 
while the coefficient Ky takes into account projection ef- 
fects, the specific mass density, the stellar orbital distribu- 
tion (such as velocity dispersion anisotropy and rotational 
support), and the effects related to the presence of dark 
matter. Equations Q and (0 imply that in real galax- 
ies, no matter how complex their structure is, T*/Ky is a 
well-defined function of any two of the three observables 
(L, i? e ,(To). For example, by eliminating <to from eqs. Q 
and (0) one obtains that along the FP 



oc R, 



(2-10/3+a)/ 'a r(5/3-a)/ 'a 



(3) 



where the dependence of the ratio T*/Ky on galaxy prop- 
erties is commonly referred as the "FP tilt" . The physical 
content of eq. © is truly remarkable: all stellar systems 
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described by eq. @ are in virial equilibrium, but only 
those for which T*/Ky scales according to eq. (0 (and 
with the same scatter) correspond to real galaxies. In other 
words, eq. |J3} indicates that structural/dynamical (Ky) 
and stellar population (Y») properties in real galaxies are 
strictly connected, possibly as a consequence of their for- 
mation process: understanding the origin of the FP tilt is 
thus of the utmost importance for the understanding of 
galaxy formation. 

A first possibility in this direction is to focus on the 
variation of a single galaxy property among the plethora 
in principle appearing in the quantity T^/Ky, while fix- 
ing all the others to some prescribed value: we call this 
approach orthogonal exploration of the parameter space. 
For instance, one can explore the possibility that a sys- 
tematic variation of Y* with L is at the origin of the FP 
tilt, while considering the galaxies as strictly homologous 
systems (i.e., with density and dynamical structures only 
differing for the physical scales, and thus Ky =const. See, 
e.g., Bender, Burstein, & Faber 1992; Renzini & Ciotti 
1993; van Albada, Bertin & Stiavelli 1995, hereafter vABS; 
Prugniel & Simien 1996). Another possibility is to enforce 
a constant Y*, and to assume that the galaxy density 
profiles, dark matter content and distribution, stellar or- 
bital distribution, and so on, vary systematically with L 
(see, e.g., Ciotti & Pellegrini 1992; Caon, Capaccioli & 
D'Onofrio 1993; Renzini & Ciotti 1993; Djorgovski 1995; 
Hjorth & Madsen 1995; Ciotti, Lanzoni & Renzini 1996, 
hereafter CLR; Graham & Colless 1997; Ciotti & Lanzoni 
1997; Prugniel & Simien 1997). 

Orthogonal explorations lead to important results, 
but, besides starting from a (more or less) well motivated 
choice of the specific parameter assumed to be responsible 
for the FP tilt, they also bring to a fine tuning problem: 
the large variation of such a parameter along the FP, nec- 
essary to reproduce the tilt, must be characterized by a 
small scatter of it at any fixed position on the FP in or- 
der to preserve the observed small thickness (e.g., Renzini 
& Ciotti 1993; CLR). Moreover, the severity of the fine 
tuning problem is strengthened by the unavoidable projec- 
tion effects associated with the three-dimensional shape of 
galaxies, if they also contribute to FP thickness. Thus, the 
interpretation of the FP cannot be limited to the study of 
its tilt only, but requires to take consistently into account 
also its thinness. 

Recently, a statistical approach to this problem, based 
on Monte-Carlo simulations and overcoming the intrinsic 
limitations of orthogonal explorations has been proposed 
(Bertin, Ciotti & Del Principe 2002, hereafter BCD). In 
this study the authors showed that, ascribing the origin 
of the FP tilt to the combined effect of luminosity de- 
pendent mass-to-light ratio and shape parameter m in 
spherically symmetric and isotropic B}l m models (Sersic 
1968), can reconcile the FP tilt with the observed large 
dispersion of m at fixed galaxy luminosity (see Figs. 5 
and 6 in CLR and Figs. 7 and 9 in BCD). Note, however, 
that in the BCD analysis the FP thickness is entirely pro- 
duced by variations from galaxy to galaxy of their physical 



properties , as a consequence of the assumption of spher- 
ical symmetry. On the other hand, elliptical galaxies are 
in general non spherical, and the quantities entering the 
FP expression do depend on the observation angle: it is 
therefore of great interest to estimate the contribution of 
projection effects to the FP thickness, and to quantify its 
physical scatter. Few analytical works have addressed this 
issue in the past (e.g., Faber et al. 1987; Saglia, Bender 
& Dressier 1993; J0rgensen, Franx, & Kjaergaard 1993; 
Prugniel & Simien 1994; JFK96; vABS), their conclusions 
pointing in the direction of a small contribution of pro- 
jection effects to the FP thickness. A different source of 
information on projection effects is also represented by the 
end-products of N-body numerical simulations (see, e.g., 
Pentericci, Ciotti, & Renzini 1995; Nipoti, Londrillo & 
Ciotti 2002ab, 2003; Gonzales & van Albada 2002). The 
impression one gets from these simulations is that pro- 
jection effects can be significant contributors to the FP 
thickness, the range spanned by the models for changing 
viewing angle being comparable to cr; n t or more. 

We explore this matter further, by extending the BCD 
approach to a class of oblate ellipsoids with non homol- 
ogous density profiles, adjustable flattening and variable 
amount of internal velocity streaming. However, in order 
to maintain the dimension of the parameter space accept- 
able we do not take into account the presence of DM ha- 
los, and the stellar mass-to-light ratio Y» is assumed to 
be constant within each galaxy. For these models all the 
relevant quantities can be expressed explicitly, thus allow- 
ing for fast numerical calculation. The paper is organized 
as follows. In Section 2 we derive the relevant properties 
of the adopted models. In Section 3 we illustrate in de- 
tail a few representative cases, focusing on the effects of 
the various model parameters on the observational prop- 
erties entering the FP relation. In Section 4 the results of 
the Monte-Carlo investigations are shown, and finally, in 
Section 5 we summarize and discuss the results. Appendix 
A collects the explicit formulas describing the model inter- 
nal dynamics, while in Appendix B we derive the expres- 
sions for the associated projected quantities. In Appendix 
C the simplest model of the family (the homogeneous el- 
lipsoid) is described in detail. 



2. The models 

2.1. 3D quantities 

In our study we use a family of oblate galaxy models with 
homeoidal density distribution, belonging to the so-called 
Ferrers ellipsoids (Ferrers 1877). The density profile is 
given by 

, s f(l-m 2 )" if0<m<l, , A s 

p(m)= P0 x^ Q . {m>1 - (4) 
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where po is the central density, n > is an integer number, 
and in cylindrical coordinates 2 (R, ip, z) the isodensity sur- 
faces are labeled by m 2 = R 2 / 'R 2 + z 2 /(q 2 R 2 ). With this 
choice R t is the model semi-major axis, while its flatten- 
ing is given by < q < 1. Note that these density profiles, 
when considered in detail, are only a rough approxima- 
tion of those of real galaxies, especially for low values of 
n. However, most of the model properties that are rele- 
vant for this study show a behavior surprisingly similar to 
that of galaxy models with more realistic density profiles 
(see Sections 3). In addition the above mentioned proper- 
ties can be explicitly written in analytic form, making the 
models suitable for Monte-Carlo simulations. 

The mass within m and the total mass of the models 
are given by 



M(m) = p R 3 t 2TTqB 3/2 , n+1 (m 2 ), 
and 



M = Po R*2ivqB 



1 



(5) 



(G) 



respectively, where B a ,&(z) = J* t a ~ x (1 — t) b ~ 1 dt is the 
incomplete Euler Beta function, B(a, b) = B Q ,{,(1) = 
T(a)T(b)/T(a + b) is the complete Euler Beta function, 
and r is the complete Euler Gamma function. 

We assume that the density profiles in eq. @ are sup- 
ported by a dynamics described by a two-integrals distri- 
bution function / = f(E, L z ) (where E and L z are the 
energy and the z component of the angular momentum of 
stars). Thus, the Jeans equations reduce to 3 



dpa 2 



R 



dz 
and 

OR 



dz ' 



P( v l 



R 



dR' 



(7) 



(8) 



where <f> is the gravitational potential, vr = v z = every- 
where, the off-diagonal elements of the velocity dispersion 
tensor vanish, and a\ = a 2 (see, e.g., Binney & Tremaine 
1987, hereafter BT). The appropriate boundary conditions 
are u\ = a 2 = on m = 1 (Ciotti 2000), and so, the for- 
mal solution of eqs. 0-10 is: 



P°R 



P-^-dz, 

dz 



where Zt = q\/ R 2 — i? 2 , and 



P(v 2 



a 



R 



= R 



( dpoj , dj>_ 
^ dR 



dR 



(9) 



(10) 



2 These coordinates are related to the natural Cartesian co- 

x/R, 



y/x 2 + y 2 , COS if ■ 



ordinate system by the relations R - 
sin p = y/R. 

3 We use symbol v for the velocity in the phase space, while 
w(x) = v is the streaming velocity as defined in eq. (B2). In 
general, a bar over a quantity means average over phase-space 
velocities. 



As it is well known, the gravitational potential of 
homeoidal systems can be obtained by evaluating a two- 
dimensional integral (see, e.g., Chandrasekhar 1969), but 
in general this integral cannot be expressed in terms of 
elementary functions. From this point of view the density 
profiles adopted here are a nice exception: their potential 
can be written explicitly (for n integer) as a finite sum of 
integer powers of R and z. Thus, from eqs. © and l|10[l 
also po~\ and p(v 2 — a R ) can be written in the same way 
(their explicit expression is given in Appendix A). 

To split v 2 into streaming motion v v = v v (that for 
simplicity we assume nowhere negative), and azimuthal 



dispersion, er 2 = (v v 
Satoh (1980) /c-decomposition: 



fc 2 (v 2 



4), 



we adopt the 



(11) 



and 



4 = °% + a - fc 2 Xv2 - 4), 



(12) 



with < k < 1. For k = no ordered motions are present, 
and the velocity dispersion tensor is maximally tangen- 
tially anisotropic, while for k = 1 the velocity dispersion 
tensor is isotropic, and the galaxy flattening is due to az- 
imuthal streaming velocity (the so called "isotropic rota- 
tor" ). In principle, by relaxing the hypothesis of a constant 
k and allowing for k — k(R, z), even more rotationally sup- 
ported models can be constructed, up to the maximum ro- 
tation case considered in Ciotti & Pellegrini (1996), where 
k(R, z) is defined so that <r 2 = everywhere 4 . 

2.2. Projected quantities 

To project the galaxy models on the plane of the sky (the 
projection plane), we employ a Cartesian coordinate sys- 
tem (xf,yf,zf), with the line of sight (los) directed along 
the zl axis, and with the xi axis coincident with the x 
axis of the natural Cartesian system introduced at the 
beginning of Section 2.1. The angle between z and zl is 
9, with < 9 < 7r/2: 9 — corresponds to the face-on 
view of the galaxy, while 9 = n/2 to the edge-on view. 
With this choice, the projection plane is (x/,y/), and the 
los direction in the natural coordinate system is given by 
n = (0, — s'm9, cos9) 5 . Accordingly, the coordinates of the 
two Cartesian systems are related by 



(13) 



and the homeoid labeled by m can be rewritten in the 
observer coordinate system as: 




m 



= VAz/ 



B 



V 



(14) 



4 Note that the important issue of the models phase-space 
consistency is beyond the tasks of this work. 

5 The los vector points toward the observer, and so positive 
velocities correspond to a blue-shift. 



4 



Lanzoni and Ciotti: Projection effects on the FP thickness 



where from now on, the symbol "~" over a coordinate will 
indicate normalization to R t , and 



A = sin 2 + cos 2 0/q 2 , 
B = (1/q 2 - l)sin0cos< 



cos 



2 



(15) 



q 2 sin 2 ( 



When integrating a model quantity along the los at given 
(xt,yt), the limits on zi are derived by setting to = 1 in 
eq. (jnj 



l-i 2 



(16) 



For examples, the surface density profile is given by 



pdzi 



poRt W) B 



g,„ + lVl-^)«+l/2, (17) 



where p(xf, yf, zi) is obtained by substitution eq. Ijf 4(1 in 
eq. (0} • The quantity £ determines the size of the elliptic 
isophotes, and their (constant) apparent ellipticity is e = 
1 — q(0). For a fixed £, the major and minor isophotal 
semiaxes are a = £ and b = q(0)£, and the associated 
circularized radius is defined by the identity ttR 2 = irab, 
i.e., Ri = y/ q{0)£- In particular, the circularized effective 
radius R c is given by 



Re = #)4, 

where £ e is the solution of the equation Mg(£ c ) 
and where the projected mass within £ is 



(18) 
M/2, 



M*(£) = / T, (lt) dxi dyi ■ 
Je'<t 



where (, ) is the standard inner product and the repeated 
index convention has been applied. The analogous quan- 
tity associated to the velocity dispersion tensor is 



c 2 = (v - v, n) 2 = aijHirij. 



(22) 



By using the two definitions above, and eqs. (B.4)-(B.5), 
the expressions for v n and a 2 are: 



-Vip cos ^ sin ( 



and 



ojt + (1 - k )(v 2 - a 2 R ) cos 2 if sin 2 0, 



(23) 



(24) 



where the last identity is obtained by using eq. I|12[) . The 
corresponding (mass-weighted) projected fields are ob- 
tained by changing coordinates in eqs. ll'-il) and (|24p. and 
then integrating on zi: 



E e (£)v p (x/, yf) e 

j: e (£)V 2 (xf,yf) 
and 

2Z e {t)al{xl,yl), 



pv n dzf, 



pvldzf, 



pa n dzf. 



(25) 
(26) 

(27) 



In general a 2 will not coincide with the velocity dispersion 
we measure in observations: in fact, in presence of a non- 
zero projected velocity field v p , the correct definition for 
this quantity is 



p ((v, n) - Up) dzi 



Y, e {l) (<J 2 p + V 2 ~v 2 p ) 



(28) 



poRt2nqB [ -,n 



1 _ (1 _ £2 )n +3/2 



We obtain 



4 = y/i - 2-l/(«+3/2) jRt = R 



(19) 



(20) 



with the effective radius of the model when seen face- 
on (or in case of spherical symmetry). As can be easily 
proved, the identity R c = \J q(0)R® is a general property 
of all axisymmetric homeoidal distributions, independently 
of their specific density profile. 

To obtain the velocity fields at (x/,y/) we integrate 
along the los their projected component on n. This is 
done by transforming the corresponding spatial velocity 
moments from cylindrical to Cartesian coordinates (see 
Appendix B). For example, the los component of the 
streaming velocity field is 



(v, n) 



Vim 



(21) 



where the last expression is derived from the identity 
(v, n) 2 — a 2 + v 2 . Note that, independently of the los 
orientation, on the isophotal minor axis yf (where, by def- 



inition, cos (p = 0) v n , Vp 2 , and v 2 vanish, and a n 



or: 

2 



on this axis a 2 is the projection of a 2 R and <r 2 os = <t_. 
In addition, the last identity holds everywhere when ob- 
serving the galaxy face-on (0 = 0), or in the case k = 0. 
Since the observed velocity dispersion is always measured 
within a given aperture, we finally integrate af os over the 
isophotes (even though <ti os in general is not constant over 
isophotes): 



M e (£)a 2 (£) = (E e a 2 os } ei 



Z e (£/)o-f os (x/,y/)dx/dy/. (29) 

£'<e 



In Appendix B we obtain the explicit expressions for a 2 , 
V 2 and their aperture values. Unfortunately, v p cannot be 
cast in algebraic form when n > 0, and so we have to resort 
to numerical integration of eq. 1)25(1 for its evaluation; the 
details are given in Section 3. 
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Fig. 1. The projected velocity fields of the (n = 3, q = 
0.3, k = 1) model, when seen edge-on (9 — n/2). 
The color contour plot shows the velocity dispersion 
field trios normalized to its maximum value trios, max — 
0.23 Rt y/Gpa. Dashed lines represent the projected rota- 
tion field Wp/cri OS;max , with numerical values labelled in 
the figure. Solid lines are the surface brightness isophotes 
[i = — 2.51og£#(^) sampled at 1 magnitude difference, 
while dotted lines are the isophotes corresponding to R c /8 
and R c . 

As a check of the exactness of the derived projected 
fields, we use a general consequence of the projected virial 
theorem (see, e.g., Ciotti 2000), i.e. 



(E fl (a 2 p + V*)) e=Rt 



= -nfWn - n 2 2 W 2 2 - ngW 3 3, 
where K^j is the kinetic energy tensor, and 



Wij = - 



px t -—d 3 x, 



(30) 



(31) 



are the components of the potential energy tensor. For our 
models, W22 = Wii, ni = 0, and the explicit expressions 
of Wn and W33 are given in Appendix A. 

We recall that a similar approach to the one pre- 
sented in this Section was adopted by vABS, who used 
a homeoidal, modified Jaffe density profile; in particular, 
they studied the projected field corresponding (in our no- 
tation) to the quantity &f os +Vp = <J^ + V£ averaged within 
0.5i? P . 



3. The model properties 

To better understand the results of the Monte-Carlo simu- 
lations presented in Section 4, here we illustrate in details 
a few representative models, focusing on the effects of the 
various parameters on the projected velocity fields, and 
on the observational quantities entering the FP. 




1os,max 



Fig. 2. The fields cti os /cti, 
model in Fig. 1 when seen at 9 
^los, max ~0.24i? t VWo- 



and v p / erios,max of the same 
tt/4. In this case 



The projected velocity dispersion and quadratic veloc- 
ity fields in eqs. (|2TI) and l(2*H|) are evaluated by using the 
explicit expressions given in Appendix B. To obtain the 
projected field v p in eq. H25J) a numerical integration on 
zi is required. For symmetry properties we restrict the 
computation to the first quadrant of the projection plane, 
that is organized with an ellipsoidal grid made of 50 uni- 
formly spaced isophotal contours. Each contour is divided 
in 50 angles, while the los length zi + — zi- (see eq. 
is divided in 100 elements. The computation of the pro- 
jected fields on this grid, by means of a double precision 
Fortran90 code, requires ~ 20 min on a 1.2 GHz work- 
station. To check the robustness and correctness of the 
code, for all the explored models (many more than those 
presented) we have verified the projected virial theorem 
given in eq. I|30|l . and we found relative errors <10 -3 . 

The illustrative cases presented here all refer to a 
model with n = 3 and q — 0.3, but their main proper- 
ties apply to the whole family of models studied in this 
paper. Due to the constancy of the mass-to-light ratio 
within each model, the mass weighted and the luminosity 
weighted quantities are coincident. 

Figure 1 shows the edge-on view of the two obser- 
vationally accessible projected fields, a\ os and v p , in the 
isotropic rotator case. As expected, the maximum value 
of (Tios is reached at the center, while <ji os is not con- 
stant on isophotes. Note that inside the ellipse correspond- 
ing to a circularized radius of R e /S (an aperture often 
used to correct <To in the FP studies; e.g., JFK96), cri os 
is constant well within 10% (in fact, better than 1%). 
As a consequence of the adopted decomposition of az- 
imuthal motions and of the edge-on view of the model, 
the projected streaming velocity field v p (dashed lines) is 



6 



Lanzoni and Ciotti: Projection effects on the FP thickness 





.05 Q.15 


0.24 0.34 0.43 0.53 0.62 

^los/^los,max 


0.72 


0.81 


0.91 U 




Fig. 3. The fields cri os /cri OS!max and -Up/erio S , m ax of the 
model in Fig. 1 when seen edge-on, in case of a sub- 
stantial reduction of the ordered motions [k — 0.5). 
The maximum value of the velocity dispersion field is 
cios.max — 0.41 R t G p$] all the comments in the caption 
of Fig. 1 apply. 



too "rigid" Satoh decomposition: however, o~\ oa is again 
nearly constant within Re/8- 



mm 


0.03 0.13 0.22 0.32 0.42 


0.51 0.61 


0.71 0.81 


0.90 1. 




os, max 








Fig. 4. The fields crios/crios.max and Vp/erios : max of the same 
model in Fig. 3 when seen at 8 = tt/4. In this case 
^los,max ~0.31i?t VWo- 



nearly vertical 6 , its value decreases towards the center of 
the galaxy, and vanishes on the yi axis. As anticipated 
in Section 2.1 the adopted density profiles, at variance 
with real galaxies, are very flat in their inner regions: 
this is clearly visible here, where the surface brightness 
\i = — 2.5 log £#(<?) drops from the center to R e by an 
amount A/i = [(2n + l)/(2n + 3)] 2.5 log 2 < 1, in contrast 
with the drop of more than 8 magnitudes for i? 1 / 4 galax- 
ies. 

In Fig. 2 we show the same model seen at 9 = tt/4: for 
obvious geometrical reasons, the lines of constant o~\ os are 
now more similar to the optical isophotes. The field cr\ os 
within R c /8 is still constant with very good approximation, 
even though R c has increased according to eq. |1 8jl. The 
field v p is deformed by projection effects, and its value, 
normalized to the maximum of trios, is lower than in the 
edge-on case, as expected. 

When the amount of ordered rotation is substantially 
reduced (for example, by assuming k = 0.5), the resulting 
velocity fields are modified as shown in Fig. 3 (9 = it/ 2) 
and in Fig. 4 (6 = it/ 4). Direct comparison with Figs. 
1 and 2 indicates that such a reduction of k moves the 
maximum of <7i os from the center to the external regions 
of the model, a consequence of the increase of o~ v at large 
galactocentric distances on the equatorial plane in order to 
sustain the model flattening. This trend of o~\ os is usually 
not observed in real galaxies, and it can be ascribed to the 

6 The n = model is an exception: its lines of constant 
v p are always parallel to the isophotal minor axis (eq. [C.7]), 
and in the edge-on, isotropic rotator case, cr\ OB is constant on 
isophotes (eqs. [C.9]-[C10]). 



While we have shown here a few illustrative examples, 
we find that for all the models studied in detail cti os stays 
nearly constant within R e /8. For example, for the n = 
model, by expanding eq. (C.ll) for £ << 1 we obtain 

^(0,0) [ 2 Ww 3 J \ 2 ,y 

and so cr p (0, 0) and <7i os ,a(-Rc/8) differ less than 0.1%, while 
for the n = 3 models represented in Figs. 1-4 the two 
quantities differ less than 1%. This implies that when us- 
ing apertures of the order of R c /8, the average los velocity 
dispersion can be safely replaced by its central value, i.e., 
oios,a(-Rc/8) — crios(0, 0) = (7p(0,0), independently of rota- 
tion and los inclination angle (note that the last identity 
holds exactly; see Section 2.2). These considerations are 
also confirmed by the cuspier models discussed in vABS. 

We now describe how models "move" in the edge-on 
view of the FP, when changing the intrinsic (i.e., n, q, and 
k) and observational (9) parameters. In Fig. [3^, we illus- 
trate the behavior of three n = models in the (log R d ,fp, 
logi? c ) space, where log i? ,FP = alogero+/3fog^e+7- Two 
models have maximum intrinsic flattening (q — 0.3), and 
differ for the amount of rotation (k = 1, and k — 0); the 
third one is a rounder (q — 0.5) isotropic rotator. Owing to 
the particularly simple expression of <7ios,a (eq. [Cll]), we 
also investigate the effect of adopting different apertures 
for the estimate of <7q. In all cases, varying the projec- 
tion angle from to tt/2, makes R c to decrease according 
to eq. I|18l) , thus producing the vertical down-shift of the 
representative models; obviously, such a shift is smaller for 
the rounder galaxy. When R e decreases I e increases, and 
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galaxy models move along straight lines of constant incli- 
nation with respect to the edge-on FP, independently of 
their specific density distribution (see comments below eq. 

and provided that oq is only weakly dependent on 
the los inclination. We find that even when o§ is measured 
within apertures of up to the order of R ei the effect on the 
model displacement in the FP is marginal, independently 
of the viewing angle and the amount of rotational support. 
The only relevant case is when the total aperture is consid- 
ered (circles) and the galaxy flattening is supported by the 
azimuthal velocity dispersion (these results can be qualita- 
tively interpreted by using eq. [C.ll] with t = l) 7 . In Fig. 
Eh>, the effect of the viewing angle for different values of 
the shape parameter n is illustrated. At variance with Fig. 
[SJl, all models have the same flattening (q = 0.3), mass, 
mass-to-light ratio, and truncation radius. The amount 
of projection effects is quantitatively similar for different 
values of n, being mainly due to the dependence of R e 
on the viewing angle. In summary, since in all cases the 
directions along which models move are not parallel to 
the FP best fit line, projection effects do contribute to the 
observed FP scatter, with effects <2<7i n t- 

Interestingly, from Fig. [SJd it is also evident that the 
trend of n along the FP is in agreement with what 
found observationally when galaxy light profiles are fit- 
ted with the i? 1/m models (e.g., Caon et al. 1993; CLR; 
Graham & Colless 1997; Prugniel & Simiens 1997; Ciotti 
& Lanzoni 1997; BCD): in fact, in this latter class, an 
increase of m corresponds to the galaxy density profile 
being radially flatter in the external regions, as Ferrers 
models behave for decreasing n. However, the amount 
of this effect in Ferrers density profiles is substantially 
smaller. This can be estimated by considering the be- 
havior of Ky — GM/i? c cTp(0, 0): in the spherical limit, 
4.93 < Ky(n) < 6 for < n < 10, to be compared with 
the range 7.96 > Ky(m) > 1.75 for 1 < m < 10 in case 
of i?, 1 /" 1 models (BCD). Thus, while structural non- 
homology alone, with constant T„, is sufficient to repro- 
duce the whole tilt of the FP in the case of the R}l m mod- 
els (CLR; Ciotti & Lanzoni 1997; BCD), this is not true 
for the Ferrers ellipsoids. Also note that, at least in the 
spherical limit, the values of Ky(n) are within the range 
spanned by the virial coefficient in the case of the R x l m 
models. In particular, they are close to the value of 4.65 
that characterizes the de Vaucouleurs profile i? 1 / 4 . Since 
Ky is the only model-dependent property that explicitly 
enters the FP relation (through eq. 0), this ensures that 
the class of models we are using is suitable for our inves- 
tigation. The behavior of Ky[n) for different flattenings 
and viewing angles is summarized in Tabled 



7 Note that the observed velocity dispersion entering the FP 
relation is usually corrected to a circular aperture with diame- 
ter 1.19 kpc (e.g., J0rgensen at al. 1999), corresponding to 
a radial range ~ 0.05 R a ~R c for h = 0.5, and for typical values 
of R c . In any case, the FP equations derived by using R a /8 
or the fixed metric aperture are in mutual good agreement 
(JFK96). 
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Fig. 5. Panel a: effects of los direction and spectroscopic 
aperture on three n = models, arbitrarily placed on the 
edge-on FP (solid line), as a function of flattening (q) and 
amount of ordered rotation (k) . Dotted lines mark cxint . In 
each model, the three aperture radii are Re/8 (squares), 
R e (triangles), whole aperture (circles), and for each aper- 
ture radius, the three points at decreasing y-axis corre- 
spond to 9 = 0, 7r/4, 7t/2. Panel b: effects of los direction 
and density profile on the position of the models along the 
edge-on FP, when <Jo = cr p (0, 0). As in panel a, 6 increases 
along the arrows. 

4. Simulations and results 

4.1. The numerical procedure 

In this Section we extend the statistical approach pre- 
sented in BCD, and we determine the most general man- 
ifold (in the parameter space) defined by the models 
that lie on the observed FP. In practice, for each seven- 
dimensional point (n, T*, L, i?t, q, k, 9) in the model space, 
we determine i? c , I c , and oo- Then, we check if the model 
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Fig. 6. Accepted models (dots) selected according to eq. 
ffity and er inp = 0.054, plotted in the edge-on (Panel a), 
and face-on (Panel b) views of the FP. For models in Panel 
a, 0mt = 0.057. Solid triangles are the Coma cluster galax- 
ies in Gunn r (from JFK96). Panel c: the edge-on FP for 
the models selected in the "zero-thickness" approxima- 
tion. Solid and dashed lines in Panels a and c represent 
eq. and its intrinsic scatter, respectively. The dotted 
line in Panel b marks the "zone of exclusion" (JFK96). 



Table 1 . The quantity for Ferrers ellipsoids as a func- 
tion of n, q, and viewing angle 9, when ctq = cr p (0,0). 



K v 




9 = 1 


q = 0.3 
9 = 


q = 0.3 

6 = tt/2 


n = 





4.93 


8.29 


15.13 


n = 


3 


5.78 


8.61 


17.36 


n = 


6 


5.92 


8.68 


17.72 


n = 


10 


6.00 


8.71 


17.89 



"belongs" to the FP. The observational FP that we take 
as a reference is the one obtained by JFK96 for the Coma 
cluster galaxies in Gunn r (eq. PQ). For its intrinsic scatter 
we adopt the value a- mt = 0.057, as quoted by JFK96 for 
the galaxies with ctq > 100 km/s. 

The domains of model parameters considered in the 
simulations are the following: 0<n<6, 1<T*<10 
(different from galaxy to galaxy, but constant within each 
model), 2.7 < L < 50 (in 1O 1O L0, the same range of 
values spanned by the JFK96 Coma cluster galaxies), 
1 < R t < 200 (in kpc), 0.3 < q < 1, and < 9 < tt/2. The 
values of n and T* are randomly extracted from uniform 
distributions, while power-law distributions p(L) oc L^ 1 
and p(Rt) oc R^ 1 ' 5 have been used to extract L and Rt 
by means of the von Neuman rejection technique. The as- 
sumption of strongly non uniform input distributions for 
L and Rt was necessary in order to end (after the FP selec- 
tion) with galaxy models having a luminosity function and 
a distribution of effective radii in agreement with those ob- 
served (see Sect. 4.3). For the extraction of the flattening 
q a fit to the observed distribution of intrinsic cllipticity 
(as derived for a population of oblate spheroids by Binney 
& de Vaucouleurs, 1981) 



0.6(q + 0.055)(<7 2 - 1.07g + 0.377) 
q 2 - 1.197^ + 0.372 1 



(33) 



has been used. Even though the assumption of a popula- 
tion made of oblate spheroids only is not fully consistent 
with observations (see, e.g., Binney & Merrifield 1998), it 
is acceptable for our investigation, and it is also consistent 
with the geometry of the adopted models. Concerning the 
estimate of the model central velocity dispersion, we recall 
that the way how models move along the FP is almost the 
same when using apertures of Re/8 or R e , and so, accord- 
ing to the results of Section 3, we assume Co = c p (0,0). 
In such a way the dimensionality of the parameter space 
is reduced by excluding k from the analysis. Finally, to 
sample the effect of the los inclination, we compute the 
projections of each model along 11 viewing angles equally 
spaced in < cos# < 1. 

For each projection angle we first check whether R e 
and Co are within the ranges 1 < i? /kpc < 20 and 
100 < <7o/(km/s) < 350; if not, the model is discarded 
as unrealistic, otherwise we construct the angle average 
quantities (logi? c ), (logero), (logic), and (logi? ,Fp) = 
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a(logCTo) + /3(logle) + 7, we calculate the quantity 8 
(AFP) = (logi? e ) — (log i? c ,Fp), and we apply the following 
criteria to check whether this candidate model "belongs" 
or not to the FP. Since the FP residuals are consistent 
with being distributed as a Gaussian (JFK96), we require 
that, according to the von Neuman rejection technique, 
the angle average model is extracted from the distribution 

p«AFP)) cx exp(-(AFP) 2 /2a i 2 np ), (34) 

where a lnp is an input parameter. We finally accept the 
model if it also belongs to the face-on FP, i.e., if it 
satisfies 0.54 a; + y<4.2, with 2.66 a; = 2.21 (log R c ) - 
0.82(log/ c ) + 1.24(logoo), and 1.49 y = 1.24(log/ c ) + 
0.82(logcr ) (JFK96). 

The end product of a complete run is the data sam- 
ple composed by the 11 projections of each accepted 
model. For this sample, in the (log R c ,fp, log Rc) space 
we estimate both the linear best-fit and the rms of the 
residuals around the best fit line, rms (AFP). The pro- 
cedure is repeated by changing the input parameters 
finp until rms (AFP) = 0.057. The physical thickness of 
the FP is then evaluated as rms ((AFP)), and the con- 
tribution of projection effects is estimated as cr pro j = 
\J rms (AFP) 2 - rms((AFP)) 2 . 

Note that in this approach, due to the high dimension- 
ality of the parameter space, to the thinness of the FP, and 
to the von Neuman rejections on L, R t , q, and (AFP), we 
usually need to calculate several hundreds thousand pro- 
jected models for each choice of CTi np . This is feasible with 
the adopted class of models, because i? e and Co can be 
expressed in a fully analytical way. 

4.2. Projection effects on the FP thickness 

Following the procedure described above, we find that for 
Cinp ~ 0.054 the sample of accepted models defines a 
synthetic FP that matches very well the observed one, 
both in the edge-on (Fig. HJl) and in the face-on (Fig. 
EJd) views. The model FP is characterized (by construc- 
tion) by rms (AFP) = a- m t, while its physical thickness is 
rms((AFP)) - 0.052 = 0.91<r int . It follows that a proj ~ 
0.41cri n t. In this simulation, the fraction of accepted mod- 
els is - 2.2%. 

Another possibility to estimate the contribution of pro- 
jection effects to the FP thickness is that of selecting 
the angle averaged models from what we call the "zero- 
thickness" FP: in practice, we adopted a lnp = 0.001 in 
eq. (34), so that the dispersion produced by the accepted 
models when seen from the 11 different los is entirely 
due to projection. Note how, in agreement with the re- 
sults shown in Fig. 5, the final data set nicely fills the la 
strip in Fig. Et- In these zero-thickness realizations, the 
fraction of accepted model is ~ 0.4%, and rms (AFP) ~ 
0.45 (Ti n t. Accordingly, we quantify the physical FP scatter 

as o-p hys = ^af nt - cr 2 roj ~ 0.89cr int . 

8 We define residuals about the FP the quantity AFP = 
log R c - log i? c , PP 



As a test of the robustness of the above estimates we 
also explored the case in which the distribution of AFP is 
a step function, i.e., instead of using eq. (|34J) . we accept 
the model if 

5 = |10< AFP > - 1| < ^p, (35) 

and we chose Si np so that rms (AFP) of the selected models 
equals a m t: the resulting values for <7 p h ys and (T pr oj are in 
perfect agreement with those obtained with the previous 
approach. 

For the whole sample of models that reproduces the 
observed FP no systematic trends with the FP residuals 
are shown by n, q, R c , L, o~o, while we find marginal cor- 
relations between the apparent ellipticity e = 1 — q{9) and 
the mass-to-light ratio T„, and AFP. In agreement with 
the analysis of JFK96 (cfr. Fig. 8a therein) and Saglia et 
al. (1993), flatter galaxies in projection are preferentially 
characterized by positive residuals (Fig. EJi), while AFP 
decreases from positive to negative values for increasing 
Y* (Fig. 03). 

4.3. The FP tilt 

We now address the issue of the FP tilt, and compare the 
properties of the models that reproduce the observed FP 
against the available observational counterparts. 

In Fig. |21 the distributions of the model properties are 
shown with solid lines, while those of the adopted obser- 
vational sample are represented with dotted lines. From 
Figs. |H1 it is apparent how the FP selection modifies the 
input power-law distributions of L and R t into distribu- 
tions that match remarkably well the observed ones for 
L and R c . In particular, the result should be contrasted 
with the one (not shown here) obtained when extracting 
L and Rt from uniform distributions: in that case, logL 
and i? c peak at ~1.1 and ~14, respectively. An interest- 
ing case is presented by the flattening q: in Fig. |HJi it is 
apparent how the input distribution (dotted curve) is not 
modified by the FP selection. Also the effect on the shape 
parameter n is not very strong, even if the FP seems to 
be marginally selective against the lowest and the high- 
est values of n (Fig. On the contrary, the effect on 
the mass-to-light ratio is remarkable: its input uniform 
distribution has been substantially altered towards small 
values by the FP selection (Fig.|Sf). Note the peak around 
T* ~ 3, a value in good agreement with the commonly ac- 
cepted stellar mass-to-light ratios in elliptical galaxies in 
the Gunn r band. 

We now describe how the model properties vary along 
the FP. While n and q do not show any particular trend or 
segregation, L and Y» are found to systematically increase 
with i? 0i FP (Fig. EJ. In both cases, the scatter is large, 
but almost constant for fixed i? e ,FP- As shown in Fig. EH, 
not only the trend of L with R c ,fp, but also the overall 
region populated by the models in this space, correspond 
remarkably well to those found in the observations. 

The distributions of L and Y* along the FP translate in 
a well defined mutual dependence of these model proper- 
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Fig. 7. Distribution of the residuals about the FP as a function of apparent ellipticity (panel a) and mass-to-light 
ratio (panel b), for models reproducing the observed FP. Points within the dashed lines correspond to models with 
|i? c /i?o,FP-l|<0.15. 



ties, as illustrated in Fig. E3 (small dots). Also in this case 
the agreement with the estimates from the observational 
data is remarkable (cfr. Fig. 3a in JFK96). The large scat- 
ter in T* at fixed L is by construction consistent with the 
small thickness of the FP: apparently, other model prop- 
erties vary within the sample of accepted models so that 
their combined effect is to maintain the FP thin. A clearer 
view of the situation can be obtained by considering only 
galaxy models lying on the idealized zero-thickness FP 
(open circles): the relation between and L is better 
defined now, even though the scatter in Y* at fixed L 
is still significant. In case of an orthogonal exploration 
based on a systematic trend of T* with galaxy luminosity, 
the set of models would be just a 1-dimensional line in 
Fig. 10 (cfr. to Fig. 6 in BCD). Here, such a case can be 
mimicked by restricting further to a sub-sample of models 
characterized by a small range of flattenings (for example, 
0.9 < q < 1, filled triangles in Fig. IHty . For this latter 
data set, the strict correlation T» oc L 3 is obtained, as 
predicted by eq. © with Ky ~ const. 




log L 



We recall that in this class of models n has only a mi- 
nor effect in determining uo (at variance with the shape 
parameter in the case of R}/ m models). However, the be- 
havior of the density profiles and L is qualitatively the 
same in the two families, as discussed at the end of Section 
3. In fact, when limiting to Ferrers ellipsoids with constant 
Y* and q, a correlation between the shape parameter and 
the luminosity appears, with n decreasing for increasing 
L. 



Fig. 10. Y* vs. L for the models that reproduce the ob- 
served FP (small dots). Circles correspond to models be- 
longing to the zero-thickness FP, while black filled trian- 
gles are their sub-sample in which 0.9 < q < 1. 

5. Discussion and Conclusions 

In this paper we have explored the importance of projec- 
tion effects on the FP thickness. We extended the sta- 



Lanzoni and Ciotti: Projection effects on the FP thickness 



11 




0.5 1 1.5 5 10 15 20 2 2.2 2.4 2.6 

logL R e logCfj 




0.4 0.6 0.8 1 1 2 3 4 5 6 2 4 6 8 10 

q n t. 



Fig. 8. Histograms of the properties of the models reproducing the observed FP (solid lines). Dotted histograms in 
panels a, b, c correspond to the observational data for the Coma cluster in Gunn r (from JFK96). The dotted line in 
panel d represents the input distribution of q, as given in eq. 



a 



1.5 



j 
o 



,5 



A A 
A A A A 



A * *A 

Aa A A 

I i f I 



: b 



0.5 1 1.5 
log R 



0.8 

-2 0.4 
0.2 




e.FP 



0.5 1 1.5 
log R e FP 



Fig. 9. Distribution along the FP of luminosity and mass-to-light ratio of the models that reproduce the observed 
FP (dots). Solid triangles correspond to the observational data for the Coma cluster in Gunn r (from JFK96). 



12 



Lanzoni and Ciotti: Projection effects on the FP thickness 



tistical approach introduced by BCD to a class of oblate 
galaxy models, with variable density profile, flattening, 
amount of rotation, and fully analytical spatial and pro- 
jected dynamical properties. In particular, after generat- 
ing galaxy models corresponding to random choices of the 
model parameters, we retained only those defining a FP 
with the same tilt and thickness as the observed one. With 
this approach not only we quantified the importance of 
projection effects on the observed FP thickness (thus also 
determining the "physical" FP scatter associated to the 
dispersion in the galaxy internal properties), but we also 
studied in a consistent way the possible origins of the FP 
tilt. Such a framework represents a valuable alternative 
to the somewhat arbitrary orthogonal explorations, where 
the property responsible for the FP tilt is selected a pri- 
ori, and a fine tuning problem for the selected parameter 
is unavoidable. In the present approach instead, we find 
that model properties can vary significantly from galaxy 
to galaxy, while preserving the observed FP thinness (a 
result in agreement with what found by BCD). Of course, 
we stress again that (at a deeper level) the very existence 
of the FP with such a small scatter is a fine tuning case, in 
the sense that stellar population (T») and structural and 
dynamical properties (K\) are tightly correlated, as de- 
scribed by eq. ©• The reason for that can only be found 
in a comprehensive theory of galaxy formation and evolu- 
tion. 

The main results of the present work can be summa- 
rized as follow: 

— For the adopted class of models the contribution of 
ordered streaming motions to the observed velocity 
dispersion is negligible when small/medium apertures 
(Ss Rc) are used for the spectroscopic observations. 
This implies that a systematic decrease of rotational 
support with increasing luminosity is not at the origin 
of the FP tilt of elliptical galaxies. However, the con- 
tribution of rotation can be non negligible when using 
larger apertures, and this must be taken properly into 
account when studying the FP at high redshift. 

— When observed from different los inclinations models 
move along a direction that is not exactly parallel to 
the best fit line of the edge-on FP, thus confirming 
that projection effects do contribute to the observed 
FP scatter. The amount of such a shift depends mainly 
on the intrinsic flattening of the models (being larger 
for more flattened systems) and, to a minor extent, on 
the adopted aperture used for the determination of <tq ■ 

— The estimated contribution of projection effects to 
the observed dispersion of logi? c around the FP is 
~ 0.4 (Ti n t, while the FP physical scatter (as determined 
by variations of the physical properties from galaxy to 
galaxy), is ~ 0.9eri nt . It follows that, when studying 
the correlations between galaxy properties required to 
reproduce the FP tilt and thickness, spherical models 
are an acceptable approximation. 

— Weak correlations of the FP residuals with galaxy 
mass-to-light ratio and apparent ellipticity are found. 



The latter is in good agreement with the results of 
JFK96 and Saglia et al. (1993). 

— For the models that reproduce the observed FP, a very 
good agreement between their properties (L, R e , <jq, 
and q) and the observational data if found. Also the 
trend of Y* with L corresponds very well to the one 
estimated from the observations. In order to get these 
results, it has been crucial to extract the models from 
steep power-law distributions of L and Rt, while uni- 
form distributions in input produce an excess of models 
with bright luminosities and large effective radii. While 
the conclusions about the contribution of projection 
effects to the physical FP thickness would remain un- 
changed, it is still unclear why the requirement that 
models belong to the FP is so selective against low L 
and R c . 

— For what concerns the FP tilt, L and T* of the ac- 
cepted models appear to systematically increase with 
Rc.fp- The corresponding increase of T* with L is 
also well defined and in very good agreement with 
the observational estimates. In addition, in the (T*, L) 
plane, the models appear to be segregated in terms of 
the flattening: at any fixed L, systems with low T* are 
preferentially rounder than those with high T*. 

— The ranges of variation of the shape parameter n, flat- 
tening, and mass-to-light ratio for models of given lu- 
minosity can be very large, although consistent (by 
construction) with a thin FP. This naturally solves 
the fine tuning problem met by the "orthogonal ex- 
ploration" approach, and it provides better agreement 
with observational data. In practice, model parameters 
mutually combine in such a way that a large dispersion 
of galaxy properties is allowed at any fixed location of 
the FP, while preserving its thinness. 

For simplicity we have used, as a guiding tool for 
the present investigation, simple one-component, two- 
integrals oblate galaxy models. The study would be best 
carried out with other families of models, better justified 
from the observational and physical point of view, but for 
the present purposes the simple models used provide an 
adequate demonstration. In fact, although the density dis- 
tributions adopted here are not a good representation of 
real elliptical galaxies, the conclusions obtained about the 
projection effects on the FP can be considered robust. In 
fact, it is also reassuring that the displacements in the 
FP space of our models due to a change in the los di- 
rection agree both qualitatively and quantitatively with 
the analogous results obtained from the end-products of 
N-body simulations. This is because the way R c (and 
thus I c ) depends on the viewing angle is the same for all 
homeoidally stratified density distributions. A substantial 
improvement of the present exploration, assuming more 
realistic density profiles, would require much more time 
expensive simulations, since numerical calculation of the 
projected velocity dispersion would be necessary. 
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Appendix A: 3D quantities 

A.l. Potential expansion and gravitational energy 
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The general quadrature formula for the potential at inter- Note that in the notation of eq. (A. 5) W\ — 4>oio(q) and 
nal points of heterogeneous homeoidal distributions, can u>3 = <fon(<z), i- e - ; 



be rewritten in the specific case of oblate Ferrers ellipsoids 
as 
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A(t) = (R^ + t) x /q*R? + T. 



(see, e.g., BT). In the limit of nearly round models (q — > 1) 
wi ~ 2(l + 2 9 /3)/5 + 0(l-q) 2 and w 3 ~ 2(3 - 4q/3)/5 + 
(A.3) 0(1 -qf. 

A.2. Velocity dispersions expansion 

(A. 4) We start by noting that eq. 10} with n integer can be 
rewritten as 



After a double expansion of identity (A.2) a cumbersome 
but trivial algebra shows that eq. (A.l) can be written as 
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where 



Pnij(q) 



(-1)' 
qZ{ % -i) \i J \j 



(A.14) 



/(?) 



q (-iy 



(n + !)(» + 1/2) 



n+l 

i 



2F1 3 



With the aid of this expansion, we can proceed to the solu- 
tion of eqs. (|9"|l- (|10JI : at first sight the simplest way could 
appear the direct substitution in these equations of ex- 
pansions (A. 5) and (A.13). On the contrary, even though 
these substitutions show that the solution of the Jeans 
equations can be obtained as finite expansions, for com- 
and a Fi is the standard hypergeometrical function defined puta tional reasons this is not the most efficient way to 

perform the calculations. For example, by using the fact 
that p and </> depend only on z 2 , in eq. @ we changed in- 



— ,i H — , i H — ; 1 — q' 
2' 2' 2' 1 



(A.6) 



by 

2 Fi(a,&, c;z) 



r(c) 



T(b)T(c-b)J 



t b t) c b zt) a dt, (Aeration variable from z to z 2 , we expanded the quantity 

dA^{r)/dz 2 in eq. (A.l), and finally we integrated the 



with c > b > and | z |< 1. Note that 2 Fi(a, 6, c; 0) = resulting expression multiplied by eq. (A.13). The result 

and that for i and 3 integers, 2F1 in eq. (A.6) reduces to is 

a combination of elementary functions, n i a k 

The components of the gravitational energy tensor pa 2 R ^GR 2 pl^^^^C nijk i(q)R 2< ' 3+l) z 2< - k ^ l \(A.15) 



used at the end of Subsection 2.2 can be evaluated from 
the useful identities 



Wi. 



GR*Tr 2 qa 2 WiS^ 



A* 2 (m)d 



(A. 



52(3+/) ?2(*-0 
i=0 3=0 k=0 1=0 

where a = 2n + 1 —i, and Cnijkl(q) = C1C2C3 with 

Ci = n q 2{n - l - k+l)+Z (1 - q 2 y , (A.16) 



where ai, 02, 0,3 are the three semi-axes of the homeoid, 
Wi are functions of the flattening q, A "J is given in eq. 



(2 



(-l)J+ fc 
2n + 1 - « 



2n + 1 - i 
k 



(A.17) 
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Ca - B ( n + - 



2F1 



i + 3,i 

3 3 

n H — , n H 

2' 2 



.(A.18) 



The r.h.s. of eq. is the sum of two contributions. The 
expansion of Rdpa^/dR is immediately obtained from eq. 
(A. 15), while following the same approach used for the 
derivation of eq. (A. 15) it can be proved that 



dtp 



2a ~2b 



Rpf^ = gr\pI E E E E (i) R za z 

' i=0 j=0 k=0 1=0 



B.2. Projection ofR 2a z 2b 

In order to integrate along the los eqs. l|2l)|) and (|T7Jl . we 

first note that costp — x/R — xr/R, i.e., the term cos 2 tp 
brings in the integrand only the quantity R~ 2 , being xi 
fixed at the chosen position in the projection plane. In 
addition, from eqs. JUJ, (A. 15), and (A. 19), it results that 
p(v 2 — crjj) is a sum of terms R 2a z? h , with a and b non 
negative integers. As a consequence, the projection of the 
quadratic velocity fields reduces to the evaluation of finite 
sums of terms of general form: 



(A.19) 



I a ,b(x', y') 



where a = j + / + 1 , b = i — j + k — I, and 
2tt {-l) l+k 



R 2a (xf, yl, zi) z 2b (xi, yl, zi) dzl, (B.6) 



Vnijkl (q) 



q 2(k-l)-i i + 3/2 



l-<t 



(A.20) 



where R and z are expressed in terms of xf, yl, and zi 
following eq. (|13fl . The result is 

a a 2b a — a c 

Iajbixi, yi) =i?tE E E E E VA & l) C(xr, yi), (B.7) 

a=0 0=0 7=0 6=0 £=0 



Note that a simple check of eqs. (A. 15) and (A.19) shows where c = a + j3 
that the cylindrical radius R 2 can always be factorized in 



the expression of the quantity p(v 2 — a 2 R ): in other words, 
the latter quantity vanishes (as expected) on the symme- 
try axis of the model. 



Appendix B: Projected quantities 

B.l. From cylindrical to Cartesian coordinates 

To project the velocity fields on the plane of the sky 
according to eqs. (|25() - (|29|l . we first transform velocities 
from cylindrical coordinates, where the Jeans equations 
are solved, to the natural cartesian coordinates, in which 
the projection is particularly simple (eqs. |21)-|22j) 



7+1, 
L2" 



c 

„2\c— 6/ 



a — a 
5 



,(B 



\2(a+j3+b)-t + l 



(cos/ 



^2(q+7+5)-e+1 

(B.9) 
(B.10) 



Vr COS tp 

vr sin tp - 



- v v sin tp, 
v v cos tp, 



y v z = v z . 
By definition 



and 



- f /(x,v,i)d 3 v, 

p J 



(Vj - Vi)(Vj ~Vj) = ViVj - ViV 



and so, following the choice of / = f(E,L z 
(B.l) 

—v v sin tp, 
v v cos tp, 
0, 



(B.l) 

(B.2) 

(B.3) 
from eq. 



^ = <7 c (l-^) c - e (sinl 
q{9)- 2 ^, 

and 

C = yi 2 ( a + b + s )-^+ 1 5/2(0-0-*) q _ giy/z^ 
with £ 2 given in eq. Q15[l. 

6.3. Integration over isophotes 

In order to derive the quantity af os a (£) we inte- 
grate T,g(£)af os over isophotes, and this requires a bi- 
dimensional integration of the quantities I a ,b{xi, yi) in eq. 
(B.6). From eq. (B.7) the only factor affected by the inte- 
gration is C, and so, with the natural parameterization: 



xl = t cos A, 
yl = q(9)t sin A, 

according to eq. I|29l) we have 

c 2k i 



(B.ll) 



(C(x/,y/))i 



C[t cos A, q(6)t sin \}q(9) t dt d\. (B.12) 



cr xx = 0% cos 2 tp ■ 
<jyy = a\ sin 2 tp - 

<Jzz = 2 R , 

r 2 „% < 



- °% sin 2 tp, 
a% cos 2 tp, 



(B.4) ^he resm t of the integration can be rewritten as 
(C) e = R 2 B(9, q) Ax (a, b, a, S, e) T x (a, b, e, I), 
where 

B=[l- (-1) £ ] g(6/)2(a+6+«+l)-e ) 



(B.5) 



V x y = (Cfl - cr 2 ) sin cos 



'xy 
<?xz 



' yz 



= 0. 



A x = B a 



a 



r 1 

5+ 2' a " 



■'-5 



i 



(B.13) 

(B.14) 
(B.15) 
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and 



B a +6- £ /2+3/2,e/2+l(^ 2 



~ 8(1 — q)/15. Equations (|25|l -l|27 () can be easily evaluated 
and, from eqs. (C.5), (C.6), Q and JHJl, it follows that 



(B.16) 



Another quantity that appears in the expression of 
°fos a W i s ( x/2 Ia,b(x!, y f ))e, due to the costp factor in eqs. 
p3")i. and jiUl. In this case, 



v p (x/,yr) 
and 



V 2 = v 2 



2ttGpqR 2 (wi — q 2 w 3 ) kx/sm9, 



(C.7) 



(C.8) 



(x/ 2 C) t = R 2 t B{9, q) A 2 (a, 6, a, 6, e) T 2 (a, 6, e, £), 
where B is given in eq. (B.14), whit 
3 



(B.17) everywhere in the projection plane, and independently of 
k. In addition, 



A 2 = B [ </ 
and 



a — 5 -\ — ,a + b + 5 



1 



B Q+ b_ e /2+5/2,£/2+l (^ 2 ) 



(B.li 



(B.19) 



al(x/,y/) = 2-KGpoRl 

V^u-? 2 ) 



+ (i - 



and, from eqs. i|28|) and (C.. 



2 \~,2 • 2 

g w 3 )xi sin i 



,(C9) 



(CIO) 



Appendix C: The n = model 



The properties of the constant density ellipsoid can obvi- 
ously be derived from the general expressions by setting 
n = 0. However, in this special case it is less cumbersome q£^2 
to derive them directly from their definition. The total 
mass is given by 



Note that from eq. (C.9) the quantity <Tp(0, 0) does not 
depend on the inclination angle 9, at variance with the 
cases n > 0: also the identities (C.8) and (C.10) are a 
peculiarity of the constant density ellipsoid. Integration 



los over the isophote i gives 



M = 



Airqpo i?3 



los, a 



iiGp R 2 A 1 + (1 - k 2 ) sin 2 Q A 2 



(C.l) 



1-(1-^ 2 ) 3 / 2 



(C.ll) 



with 



while the surface density profile in the projection plane ^ _ 2q 2 w fl — (1 — £ 2 ) 5 / 2 ] 
(x/, yi) is given by 



2p Rtq 



and the projected mass within isophote I is 
PoRl^q , 



M e (£) = 



.[!_(!_ ^3/2], 



(C.2) 



(C.3) 



From eq. lO, i?° ~ 0.608i? t - The potential inside the 
ellipsoid is 



<p(R, z) = —irGpoR 2 (wo — w\R 2 — ws'z 2 ), 



(C.4) 



and 

A 2 = ( Wl -q 2 w 3 )[2- (2 + ?-3£ 4 )Vl-^ 2 ]. 
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